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ABSTRACT 

We investigate the stellar populations of a sample of 162 Lya emitting galaxies (LAEs) at z = 3.1 in the 
Extended Chandra Deep Field South, using deep Spitzer IRAC data available from the GOODS and SIMPLE 
surveys to derive reliable stellar population estimates. We divide the LAEs according to their rest-frame near- 
IR luminosities into IRAC-detected and IRAC-undetected samples. About 70% of the LAEs are undetected 
in 3.6 fim down to = 25.2 AB. Stacking analysis reveals that the average stellar population of the IRAC- 
undetected sample has an age of ~ 200 Myr and a mass of ~ 3 x 10 8 M©, consistent with the expectation that 
LAEs are mostly young and low-mass galaxies. On the other hand, the IRAC-detected LAEs are on average 
significantly older and more massive, with an average age > 1 Gyr and mass ~ 10 10 M©. Comparing the IRAC 
colors and magnitudes of the LAEs to z ~ 3 Lyman break galaxies (LBGs) shows that the IRAC-detected LAEs 
lie at the faint blue end of the LBG color-magnitude distribution, suggesting that IRAC-detected LAEs may be 
the low mass extension of the LBG population. We also present tentative evidence for a small fraction (<~ 5%) 
of obscured AGN within the LAE sample. Our results suggest that LAEs posses a wide range of ages and 
masses. Additionally, the presence of evolved stellar populations inside LAEs suggests that the Lya luminous 
phase of galaxies may either be a long-lasting or recurring phenomenon. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: high-redshift — galaxies: 
stellar content 



1. INTRODUCTION 

Lyman-alpha emitting galaxies (LAEs) are important trac- 
ers of galaxy formation. The Lya emission is produced by 
on-going star formation in the galaxies, and the line emission 
enables discovery of objects that may be too faint to be seen in 
the continuum. LAEs therefore offer an opportunity to probe 
the faint end of galaxy formation at high redshift, and may 
serve as building blocks of larger galaxies in a hierarchical 
universe. 

In recent years, there have been great advancements in 
narrow-band surveys of LAEs. Large samples of candi- 
date L AEs now exist at redshifts z = 3-7 (e.g. Ih"u et all 
2004t iMalhotra & Rhoadsl 12004 iTaniguchi et all 120051: 
Shimasaku et al.ll2006t iDawson et al.ll2007t iMurav am a et all 
20071: lOuchi et alj|2007h . One outstanding question regarding 



the LAEs is their stellar population, which is important for un- 
derstanding the physical nature of the LAEs and their connec- 
tion with other high redshift galaxy populations such as Ly- 
man break galaxies (LBGs; IShaplev et al. 2003; Ando et alj 
[20061; iPentericci et alj 120071 IStanwav etaHl2007l) . A diffi- 
culty often encountered when studying the stellar population 
of high-redshift galaxies is the necessity to have observations 
at or redward of rest-frame optical in order to constrain the to- 
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tal stellar mass. At z > 3, this means having either deep near- 
IR or Spitzer IRAC observations. Much progress has been 
made recently, however, and several studies have suggested 
that LAEs are small galaxies with masses < 10 9 ~ 10 10 M© 
dGawiser et al.ll2006atlLai et al.ll2007l;lFinkelstein et alj|2007l; 
iPirzkal et alJl2007tlNilsson et alj2007l) . 

In this study, we focus on a sample of 162 LAEs at z = 3. 1 
discovered in a narrow-band sur vey of the Extended Ch andra 
Deep Field South (ECDF-S) bv lGronwall et all d2007l) . The 
narrow-band 4990 A imaging covers an area of 992 arcmin 2 
and reaches a completeness limit of ~ 1.5 x 10~ 17 erg cm" 2 
s" 1 (equivalent to a Lya luminosity Ll vq = 1 .3 x 10 42 erg s -1 ). 
Spectroscopic results confirm the robustness of the sample 
(detai ls of the spectroscopy will be presented in Gawis er" st all 
l2007l and Lira et al., in preparation). Of the 52 candidates with 
sufficient S/N in the spectra to yield redshifts, 5 1 were con- 
firmed to be z = 3.1 LAEs (the remaining object is a z = 1.60 
AGN with C III] 1909 A emission). The most likely contam- 
inants in our sample are z = 0.34 [O II] emitters. However, 
there should be a negligible number of these objects in our 
sample because by requiring the observer's frame equivalent 
width to be > 80 A, we eliminate all but t he strongest (and 
rarest) [O II] emitters dGronwall et al.ll2007l) . 

The narrow-band imaging in ECDF-S is supplemented 
by optical and nea r-IR observ ations from MUSYC 
dGawiser et al l l2006bl) GOODS dDickinson et al.1 120031) . 
and GEMS dRixetalJ [2004) surveys, and Spitzer obser- 
vations from GOODS and SIMPLE 8 (Damen et al., in 
preparation). Using the multi-wavelength data available in 
ECDF-S, we study the rest-frame UV to near-IR properties 
of this large sample of LAEs. We place special emphasis 
on the new Spitzer IRAC observations, which sample the 
rest-frame 0.9-2 fim emission from the z = 3.1 LAEs and 
provide valuable constraints on their stellar population and 
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FIG. 1. — R — m^ f, vs. R color-magnitude diagram of the IRAC- 
uncontaminated LAE sample. The dotted line represents the 1113 (, = 25.2 flux 
cut used to define the IRAC-undetected sample (shown as upper limits). Ob- 
jects that are undetected in both R and 3.6 (im (15 in total) are omitted. Note 
also that there are two objects that are detected in 3.6 fim but not in R (shown 
as lower limits near the right of the plot). Red R-m^ 6 > 2 objects are marked 
by blue squares, with the x-ray detected AGN candidate shown with a filled 
square. The red five-point stars are results of stacking (in order of increasing 
R) the IRAC-detected and IRAC-undetected samples. 

star formation history. 

We assume a ACDM cosmology with f2 M = 0.3, fl\ = 0.7, 
and h = 0.7. All magnitudes are in the AB system. 

2. IRAC OBSERVATIONS OF z = 3. 1 LAES 

When working with deep IRAC observations, source con- 
fusion is often an issue because of the broad IRAC PSF (~ 2" 
FWHM). Since the LAEs are in general quite faint in the 
IRAC bands, it is important to select a subsample whose 
IRAC photometry is not significantly contaminated by nearby 
objects. Beginning with the original sample of 162 LAEs, 
we select those with no detected neighbors within a 4" ra- 
dius (~ 2x the IRAC FWHM), and additionally reject LAEs 
with very bright neighbors at > 4" radius by visual inspec- 
tion. This results in a high confidence IRAC-uncontaminated 
sample of 76 LAEs. 

We further separate the IRAC-uncontaminated sample 
into an IRAC-detected sample and a corresponding IRAC- 
undetected sample. This division separates the rest-frame 
near-IR luminous objects from the less luminous ones. Of 
the 76 IRAC-uncontaminated LAEs, 52 are undetected in 
the IRAC 3.6 fim channel down to a flux limit of 0.3 /iJy 
(TO36 = 25.2). This is nominally the 2a detection limit of 
the survey. However, since the LAE positions are well de- 
termined from the narrow band imaging, it is possible to do 
slightly better with aperture photometry. Aperture corrected 
fluxes are measured from the IRAC images using 2 "5 diam- 
eter apertures. The flux cut of 0.3 /iJy corresponds approxi- 
mately to a 3er detection limit for the individual objects. The 
initial IRAC-detected sample consists of 24 objects, but 6 ob- 
jects are rejected as possible AGN or dusty objects (see be- 

8 http://www.astro.yale.edu/dokkum/SIMPLE/ 





1 1 1 1 1 1 1 1 




| | K - 7W, 6 > Z L AcS 




* AGN/QSOs 




• 8/(mLBGs 




* IRAC-detected LAEs ; 


X 




* 




; ■ J 






; 


1 SK 

: . t _ 










1 , , , 1 



18 



20 



22 
m 



24 



26 



FIG. 2. — »J3.6 — "J8.0 vs. mg.o color-magnitude diagram of the R— TO3.6 > 2 
LAEs (blue squares; the x-ray detected object is marked with a filled square), 
Z ~ 3 AGN/QSOs (green stars; Steidel et al. 200 1 Magdis et al., in p repa- 
ration), and 8 fim detected LBGs (black dots; Rigopoulou et al. 2006). The 
limits obtained by stacking the IRAC-detected sample are marked with the 
red star. 



low), leaving a total of 18 objects in the IRAC-detected sam- 
ple. 

IRAC samples rest-frame 0.9-2 fim at z ~ 3 and provides 
important constraints on the mature stellar population in the 
LAEs. In particular, the R—mse color, shown for the IRAC- 
uncontaminated sample in Figure Q] probes the relative con- 
tribution from young and mature stellar populations, and is 
therefore an age indicator. The dotted line represents the 
mi .6 = 25.2 flux limit for the IRAC-undetected sample. All ob- 
jects in this sample will fall as upper limits on the dotted line. 
The survey limits imply that we are not sensitive to young stel- 
lar populations with R-m-i^ < 0. However, stacking analysis 
can be used to obtain constraints on the average properties of 
the IRAC-undetected sample (see next section). The IRAC- 
detected LAEs typically have red rest-frame UV to near-IR 
colors, < R — W3.6 < 2, with a median of R-m^ b = 1.2. 
The red R- 1113 6 colors of the IRAC-detected LAEs imply the 
possible presence of a significant mature stellar population in 
LAEs, a point we will examine more closely in the following 
section. 

While the R-mjd color traces the age of the stellar popu- 
lation, the 3.6 /im magnitude probes the mature stellar mass 
in the LAEs. Therefore, the IRAC-detected sample, which 
has 3.6 /im magnitudes ranging from ~ 23.5-25, should 
in general be more massive than the m^.e > 25.2 IRAC- 
undetected sample. The IRAC-detected sample is also in gen- 
eral brighter in R, with the majority of IRAC-detected LAEs 
being brighter than R ~ 26. 

We also explore the possibility of a correlation between 
Lya luminosity and IRAC flux. The mean Lya luminosity 
of the IRAC-detected sample (log(LLya) = 42.47 erg s" 1 ) is 
about 40% higher than that of the IRAC-undetected sample 
(log(LLya) = 42.34 erg s -1 ). This is only a weak correlation, 
considering that the average 3.6 fim flux of the IRAC-detected 
sample is 8 x higher than the IRAC-undetected sample (Fig- 
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ure[3]l. There is also significant overlap in the Li, ya distribu- 
tions of the two samples, in part due to the large (~ 0.2 dex) 
dispersion in Li^ a within the subsamples. 

It is interesting to note that there are four objects in the 
sample, correspondin g to LAE 1, 11, 117, and 144 in the 
iGronwall et alJ (120071) catalog, with very red colors R—ms.f, > 
2. LAE 1, the brightest Lya source in the sample, is detecte d 
in x-ray and is presumably an AGN (IGronwall et al.l 120071) . 
The similar R—m^ colors of the remaining three x-ray unde- 
tected objects suggest that they may possibly be AGN as well. 
LAE 1, 11, and 117 are spectroscopically confirmed to be at 
z ~ 3. The x-ray detected object (LAE 1) shows C IV emis- 
sion in its spectrum, while LAE 11 and 117 show only Lya 
emission (Lira et al., in preparation). LAE 144 is observed, 
but no line was detected in its spectrum. LAE 1 1 and 1 17 are 
in the central part of CDF-S where the 0.5-8 keV (rest-frame 
0.1 - 2 keV) flux limit of < 3 x 10~ 16 erg s" 1 cm" 2 implies an 
upper limit of < 2.5 x 10 43 erg s" 1 for the objects' AGN x-ray 
luminosities. The four objects have R — m.4.5 colors ranging 
from 2.4 < R- m.A.% < 4, consistent with the colors of mid-IR 
selected AGN dHickox et alj|2007l) . However, the observed 
R - TO4.5 colors can also be produced by a dusty stellar con- 
tinuum with E(B-V) ~ 0.5. These four objects are the only 
8.0 fim detected sources in the IRAC-uncontaminated sam- 
ple and with m%a ranging from 21-22 they are more than 
2 mags brighter than the average IRAC-detected LAE (Fig- 
ure|5]l. Two of the four objects have m^-man co lors consis- 
tent with z ~ 3 AGN and QSOs (Steidel et alj2003b Magdis et 
al., in preparation). However, Figure|2]shows that the four ob- 
jects have 8.0 /i m magnitudes more similar to that of the 8 ^im 
detected LBGs dHuang et all 120051: iRigopoulou et al.l 2006). 
There are a total of two x-ray detected objects in the original 
sample of 162 LAEs (the other x-ray detected object is not 
in the IRAC-uncontaminated sample), giving an AGN frac- 
tion of ~ 1%. The evidence presented here cannot rule out 
the possibility that the other three red R-m^^ LAEs are dusty 
LBG-like objects, but if they also turn out to be AGN, then 
the AGN fraction is ~ 5%. 

Since the four R - m^s > 2 LAEs have rather different 
properties than the rest of the LAE population, they are ex- 
cluded from the conservative IRAC-detected sample. Two 
more objects are rejected since they are detected in 3.6 /im 
but not in R, which may imply a large amount of dust red- 
dening inconsistent with the rest of the IRAC-detected popu- 
lation. The subsequent analysis will be focused on the con- 
servative IRAC-detected sample of 18 LAEs. However, we 
will also present results derived using the larger sample of 23 
IRAC-detected LAEs (the x-ray detected LAE AGN is still 
excluded). 

3. STELLAR POPULATION 

Since the majority of the LAEs in our sample are unde- 
tected in the IRAC images, we perform a stacking analysis in 
order to measure the average properties of the sample. Again, 
we work only with the sample of 76 IRAC-uncontaminated 
LAEs. This is especially important when stacking the IRAC- 
undetected subsample because any flux contributions from 
bright neighbors could significantly skew the average. Be- 
cause of crowdedness in the IRAC images, it is also necessary 
to clean the images before stacking by subtracti ng detected 
sourc es using a PSF fitting routine (StarFinder; Diola itTet al.l 
2000). The optical and near-IR images are stacked in the same 
manner as the IRAC images, except that the subtraction of 
detected sources is not necessary in this case. Uncertainties 
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FIG. 3 . — Stacked SEDs and best-fit models of the IRAC-detected (squares, 
red line) and IRAC-undetected (diamonds, blue line) samples. The models 
assume constant SFR and solar metallicity. Only the bands with detections in 
the stacked images (S /N > 3) are included in the plot. 



in the stacked photometry are estimated by bootstrapping the 
original sample, and then repeating the stacking and flux mea- 
surement on the bootstrap samples. The uncertainties thus 
derived take into account both random fluctuations and the 
intrinsic l uminosity dispersion amo ng the LAE population. 

We fit lBruzual & Chariot] (120031) stellar population synthe- 
sis models with a Salpeter initial mass function to the stacked 
photometry. The V band flux is corrected for the contribu- 
tion of the Lya emission before the fits. The best-fit model 
to the IRAC-undetected sample is shown together with the 
stacked SED in Figure [3] Assuming a constant SFR and so- 
lar metallicity, the IRAC-undetected sample is best fitted by a 
model with E(B-V) = 0.0^, age = 160+J|g Myr, and mass = 
3^xlO 8 M . 

The average mass of the IRAC-undetected LAEs is smaller 
than that of typical LBGs by more than a factor of 30. The 
average SFR derived from the best-fit model is 2 M Q yr" 1 , 
which is not very large when compared to other populations 
of galaxies at z ~ 3 such as the distant red galaxies (DRGs; 
iFranx et al.ll2003l) which typic ally have SFR > 100 M yr -1 
(Forster Schrei ber et al.1 120041) . However, the LAEs have a 
specific SFR (defined as SFR per unit mass) of 7 x 10~ 9 yr" 1 , 
which is among the highest observed of galaxies at z ~ 3. 
In comparison, LB Gs have specific SFR ~ 3 x 10~ 9 yr" 1 
( Shanl ev et al.l2001l) and DRGs ha ve specific SFR - 1 x 10" 9 
yr" 1 (|Forster Schrei ber et al.l l2004). This means that the stel- 
lar populations in LAEs are increasing in mass at a much 
faster rate than other populations of galaxies at similar red- 
shifts, suggesting that they are in the early stages of forma- 
tion. 

The IRAC-detected sample appears to be significantly dif- 
ferent from the IRAC-undetected sample (Figure |3). Again 
assuming a constant SFR and solar metallicity, we find that 
the average stellar population of the IRAC-detected LAEs 
is best fitted by a model with E(B - V) = O.O^q, age = 
1.6 ± 0.4 Gyr, and mass = 9 ± 3 x 10 9 M . Compared to 
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the IRAC-undetected sample, the SFR of the IRAC-detected 
LAEs is slightly higher at 6 Mq yr" 1 . The specific SFR, how- 
ever, is about 10 x lower at 7 x 10~ 10 y r" 1 , similar to LBG s 
with comparable stellar population ages (Sha plev et alj2001l) . 
We emphasize that although the IRAC-detected and IRAC- 
undetected samples have distinct properties on average, the 
individual LAEs ex hibit a continuum of p roperties, as evident 
in Figure[T](see also lGronwall et al.l2007l) . If we include the 5 
red R—m-i.6 objects originally excluded from our conservative 
IRAC-detected sample (the x-ray detected LAE AGN is still 
excluded in this analysis), the best-fit age and mass increase 
slightly to 2.1 Gyr and 1.8 x 10 10 M , respectively. 

The best-fit parameters presented above do not depend sen- 
sitively on the assumed metallicity. Fitting the data with a 0.2 
Zq model increases the best-fit mass and the age of the IRAC- 
undetected sample by 30% and 60%, respectively, while the 
effect on the IRAC-detected sample is even smaller at less 
than 10%. The best-fit parameters also do not depend sensi- 
tively on the assumed star formation history. Comparing con- 
stant and exponentially declining SFR models with the same 
metallicity, we find that the two models in general give best- 
fit masses and ages that are within ~ 20% of each other. This 
is because the best-fit exponential SFR models have mostly 
long star formation timescales (> 1 Gyr) and are therefore 
very similar to a constant SFR model. 

At z = 3 . 1 , AGB stars may contribute significantly to the 
flux in the IRAC bands, which correspond to rest-frame 
0.9 - 2 /im. In order to examine how the effects of AGB 
stars wo uld change our c onclusions, we perform alternate 
fits using iMaras ton (2005) stellar population synthesis mod- 
els, which include larger contributions from AGB stars than 
iBruzual & Charlol ([2003) models. Fitting exponentially de- 
clining SFR models with solar metallicity, we find for the 
IRAC-undetected sample a best-fit model with age =100 Myr, 
star formation timescale = 250 Myr, E(B - V) = 0, and mass 
= 3 x 10 8 M Q . For the IRAC-detected sample, the best-fit 
model has age = 700 Myr, star formation timescale = 2 Gyr, 
E(B-V) = 0, and mass = 5 x 10 9 M Q . The Maraston mod- 
els in general give lower mass and age estimates, but despite 
these differences, the distinction between the IRAC-detected 
and IRAC-undetected populations remains significant. 

4. DISCUSSION 

There are several previous studies on the stellar popula- 
tion of LAEs at z ~ 3 and bevond dGawiser et al.j |2006a: 
Finkelstein et all 120071; iPirzkal et al.1 120071: iNilsson et alj 
20071) . These previous studies found that LAEs have masses 
ranging from 10 6 - 10 9 M , and ages ranging from 0.005- 1 
Gyr. Our best-fit average age of 160 Myr and mass of 
3 x 10 s Mq for the IRAC-undetected sample are broadly 
consistent with these previous studies. In a further analy- 
sis of the IRAC-und etected sample presented in this paper, 
iGawiser et alj d2007l submitted) fit two-component models 
and find that the IRAC-undetected LAEs have a total mass of 
1 x 10 9 Mq, with a young stellar component (age ~ 20 Myr) 
accounting for ~ 20% of the total mass. These results are con- 
sistent with ours which are derived using single-component 
models and should be interpreted as the average properties of 
the entire stellar population within the LAEs. 

The IRAC-detected LAEs are more luminous than the 
IRAC-undetected LAEs in both rest-frame UV and near-IR, 
and represent the massive end of the LAE mass spectrum. 
Potentially, the IRAC-detected LAEs may provide a link be- 
tween the LAEs and other galaxy populations. For instance, 
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FIG. 4. — A comparison of the R-mj f, colors and 3.6 /im magnitudes 
between 8 /.tm detected LBGs (green square), regular LBGs (black dots), and 
the IRAC-detected LAEs (red stars). The solid line shows the R < 25.5 se- 
lection limit for LBGs. The data for the regular LBGs come from IRAC 
observations of the Steidel et al. 12003) LBG sample (Magdis et al., in prepa- 
ration), and t he 8 )im LBG data co me from combining the samples of Magdis 
et al. and Rigopoulou et al. 1 2006). The mean colors and magnitudes of each 
sample are plotted as large filled symbols, with error bars showing the dis- 
persion within the samples. We additionally plot the stacked photometry of 
the IRAC-undetected sample, shown as the red filled red star with the faintest 
3.6 fj,m magnitude. 

they may be the z ~ 3 analogues to t he similarly m assive 
IRAC-detected LAEs found at z ~ 5.7 (lLai et al.ll20"07h . In 
Figure 2J we compare the R—m.3.6 colors and 3.6 /im mag- 
nitudes of the LAEs to 8 /im detected LBGs and regular 
LBGs at z ~ 3. The R — m ^ color and 3.6 urn m agni- 
tude correlate with stellar mass (Ri gopoulou et al.l2006l) . with 
the 8 /xm LBGs being the most massive (> 10 11 Mq), fol- 
lowed by LBGs (~ 10 10 - 10 11 Mq), IRAC-detected LAEs 
(~ 10 10 M Q ), and finally IRAC-undetected LAEs (~ 10 8 
Mq). The IRAC-detected LAEs occupy the faint and blue 
end of the LBG color-magnitude distribution, suggesting that 
they may be the lower mass extension of the LBG popula- 
tion. This interpretation is supported by the fact that the in- 
ferred mass of the IRAC-detected LAEs (~ 10 10 Mq) is at 
the l ow end of the mass range found for LBGs (~ 10 10 - 10" 
MQ; lPapovich et alJl200UIShaplev et al]|2001l) . Furthermore, 
most LAEs in the present sample have (observer frame) op- 
tical colors that would allow them to be selected as LBGs, 
although the LAEs tend to be fainter in th e rest-frame UV 
dGronwall et alj|2007t IGawiser et alj|2006al) . It has also been 
observed that ~ 20% - 25% of LBGs at z ~ 3 exhibit Lya 
emission strong enough to be narrow band excess objects 
dSteidel et alj|2000t IShaplev et alj|2003l) . There is thus accu- 
mulating evidence suggesting that the IRAC-detected LAEs 
may be the bridge connecting the LAE and LBG populations. 

Using the stellar mass estimates derived in the previous 
section and assuming a survey volume of 1.1 x 10 s Mpc 3 
dGronwall et alJl2007t) . we find stellar mass densities of 0.3 ± 
0.3 x 10 6 and 3 ± 1 x 10 6 M Q Mpc" 3 for the IRAC-undetected 
and IRAC-detected populations, respectively. The errors in 
the stellar mass densities include uncertainties in the mass 
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from the fits and uncertainties in the number density com- 
ing from Poisson fluctuations and a ~ 20% cosmic variance 
given the observed LAE bias of ~ 1.7 (tGawiser et al.ll2007l 
submitted). We should stress that the above stellar mass den- 
sities only account for LAEs with Lya luminosities above our 
survey completeness limit of 1.3 x 10 42 erg s" 1 . The IRAC- 
undetected LAEs account for about 9 ± 10% of the total stellar 
mass in LAEs, even though they make up 2/3 of the total by 
number. Compared to LBGs and DRGs at z ~ 3, which have 
stellar mass densities of ~ 1 x 1 7 and 7 x 10 6 M Q Mpc" 3 
respectively (Grazian et al. 2007), the stellar mass contained 
in LAEs is smaller, but not insignificant. However, it is im- 
portant to keep in mind that there may be substantial overlap 
between the LAEs and LBGs. A better understanding of the 
overlap between these two populations is necessary before a 
direct comparison of the stellar mass densities can be made. 

The present results show that the LAEs posses a wide range 
of masses and ages, from the massive and evolved IRAC- 
detected LAEs to the young and small IRAC-undetected 
LAEs. The range of photometric properties shown in the LAE 
sample suggests that LAEs exhibit a continuum of properties 
between these two extremes. The presence of both young and 
evolved stellar populations within the overall LAE popula- 
tion implies that the Lya luminous phase of galaxies may last 
> 1 Gyr, or that the Lya lum inous phase is recurring. Inter- 
estingly, IShapley et all (1200 ll) found that LBGs with best-fit 
stellar population ages > 1 Gyr also show strong Lya emis- 



sion in their spectra. The authors suggest that vigorous past 
star formation has destroyed and/or expelled the dust inside 
the galaxies, allowing the Lya photons to escape. One char- 
acteristic of these evolved Lya emitting LBGs is that they 
tend to have more quiescent SFR than the rest of the pop- 
ulation. The IRAC-detected LAEs, with their evolved stel- 
lar population and low specific SFR, are consistent with this 
scenario. If the LAE phase is recurring, then the young age 
of the IRAC-undetected LAEs implies that their number and 
stellar mass densities can be as much as a factor ~ 10 higher 
(very roughly the ratio of the age of the universe to the stellar 
age). Similarly young stellar popula tions (< 100 Myr) have 
also been found a t higher redsh ifts (Finkel stein et alJl2007t 
iPirzkal et aLll2007b IVerma et al.ll2007l) . The stochastic nature 
of galaxies with short lifetimes suggests that there may be a 
related population of undetected pre or post-starburst galaxies 
that may contribute significantly to the stellar mass and star 
formation rate densities at z > 3. 
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